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We set up simple harmonic lattice models for elastic fluctuations in bcc and fee lattices and 
the excitation of dislocations and disclinations. From these we derive, in a lowest approximation, 
universal formulas which predict melting temperatures in good agreement with the experiments. 
This new theory is more precise than Lindemann's rule by factor 2, and more predictive, since the 
size of the Lindemann number has to be fixed by experiments. In addition, our theory allows for 
systematic improvements. 
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Melting transitions 0, 0, 0, 3 are important phenom- 
ena of both technical and theoretical interest [H 
H EH El El El A hundred years ago, Sutherland 
E3 found an empirical rule, that the product bT m Mi is 
nearly constant for metals, where T m is the melting tem- 
perature, b the mean coefficient of expansion, and M the 
atomic mass. In the following year, he advanced a kinetic 
theory of solids 15], in which he postulated that melting 
would occur when the space between atoms reaches a cer- 
tain value relative to the atomic diameter, so that atoms 
are able to escape from the imprisonment by their neigh- 
bors. He carried out a hard-sphere simulation experiment 
by shaking a box containing one layer of marbles, and 
noting by removing marbles one at a time that the ap- 
parent solid-liquid transition occurs when a free volume 
of 25-33% is reached. From the vibration amplitude and 
the thermal kinetic energy, together with the above em- 
pirical melting rule, Sutherland |15| calculated the period 
of vibrations of metal atoms at the melting temperature, 
and found that the ratio of vibration amplitude to atomic 
spacing is nearly the same for all elements at melting. 

In 1910, Lindemann E3 combined Sutherland's ideas 
with the Debye theory of specific heat in solids, and de- 
rived his famous rule according to which a quasi-universal 
parameter, the Lindemann number, should be the same 
for all melting transitions. The number is usually stated 
in the form 

L = Bv liz {MlT m ) 1 l 2 

where 9 is the Debye temperature and v the volume per 
atom. Subsequent extensive tests of this rule were car- 
ried out by Gschneidner E3> Ross E3j and Crawford 
|l9| who found, however, that Lindemann's rule is not 
very reliable. First, the Debye temperature cannot be 
precisely specified since the Debye theory of specific heat 
is itself approximate. Second, extracting 8 from the data 
by a best fit of the Debye theory, the Lindemann number 
varies considerably from material to material, as empha- 
sized by Wallace |2ij . 

The most unsatisfactory feature of the Lindemann rule 
is that its size is unknown and has to be extracted 



from an average of many melting transitions. Thus, 
the Sutherland-Lindemann approach to melting has re- 
mained a purely phenomenological rough description of 
the transition. In spite of its roughness, it is often used 
even today in many works to estimate melting tempera- 
ture, this being due to a lack of a more precise melting 
theory. 

During the past two decades, several other melting the- 
ories have been developed. One is based on density func- 
tional [2ll |22j, and has provided us with an alterna- 
tive insight into the melting transition. This work, how- 
ever, describes satisfactorily only the long-wave length 
limit of fluctuations, and neglects the important role of 
the lattice structure, thus requiring essential corrections, 
which sometimes have been inserted these heuristically 
[23L l24j . 25] . Another theory of melting is based on a Lan- 
dau expansion in a symmetric tensor order parameter of 
rank four |26| . This theory is purely phenomenological 
and thus only descriptive, not predictive. 

For a proper description for the melting transition, it 
seems necessary to incorporate information of the lat- 
tice structure and the associated crystal defects into the 
theory. About 15 years ago, one of the authors (HK) 
did this successfully by constructing lattice models of 
the melting transition 0, 0. The important progress 
of these models [2?j was that they started out from the 
lattice version of the elastic lattice energy which allow for 
the thermal creation and annihilation of dislocations and 
disclinations and keeps track of their lattice elastic en- 
ergy |27l ] . They are included by means of discrete- valued 
defect gauge fields, the disclinations being essential for 
explaining the first-order nature of transition [28| ] , which 
otherwise would be of second order, as in the vortex-line 
induced A-transition of superfluid helium. 

In this respect, the older models went beyond those of 
more recent authors 0> HH wn0 considered only the pro- 
liferation of dislocation lines which, if properly treated, 
would have produced only second-order phase transitions 
pfflj ]. thus being unable to describe the melting transition. 

A simple universal melting formula was obtained from 
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the above lattice defect models in a lowest-order ap- 
proximation which combines the properties of the high- 
temperature expansion of the defect contributions with 
the low-temperature expansion of the free energy density. 
From the intersection of these two curves, one can obtain 
a formula for the melting temperature in accordance with 
Lindcmann's rule, but with a prediction of the absolute 
size of the Lindemann number Q • 

This successful theory has, however, one important 
drawback: For simplicity, it was constructed only for 
the physically somewhat exotic case of simple cubic lat- 
tices. This makes it strictly applicable only to very few 
materials — there are almost no simple cubic lattices in 
nature. If the derived formulas are applied to other cu- 
bic crystal structures such as body center cubic (bec) 
lattices and face center cubic (fee) lattices, the accuracy 
must necessarily be bad. Hence there is an obvious need 
to generalize the theory to physically more prolific crystal 
structures. 

The purpose of this letter is to present such a gener- 
alization of the melting model to fee and bee lattices. 
When generalizing the previous melting model to other 
cubic lattices, We shall see in the high-temperature limit, 
that for all cubic crystals the free energy density is only 
related to the elastic constants of the crystals and does 
not depend on the lattice structure of the system. This 
is in contrast to the low-temperature limit, where the 
free energy density involves integrals over the lattice mo- 
menta Ki and Ki which strongly depend on the detailed 
lattice structure of the crystals. It is therefore important 
to give a proper description of the low-temperature limit 
of the free energy density for different lattice structures. 

The lattice model for the cubic crystals will be based 
on the elastic energy 
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where n is the number of atoms per unit cell, em are 
oriented link vectors of cubic lattice pointing from any 
lattice site to its nearest neighbors along the positive di- 
rection, and Ui denote the Cartesian displacement field, 



For fee lattices, there are six link vectors: 

e (1) = (1,1,0), e (2) = (1,0,1), e (3) = (0,1,1) 

e (4) - (1,-1,0), e (5) = (0,1,-1), e (6) = (-1,0,1), (2) 

each lattice site x having twelve nearest neighbors at x ± 
ema/2. For bec lattices, there are four link vectors 



= (1,1,1), e (2 ) = (1,1,-1), 
e(3) = (1, -1, 1), e (4) = (-1, 1, 1), 



(3) 



and each lattice site possesses eight nearest neighbors. 
The symbols Vm and V(i) denote lattice derivative, 
defined by Vmf(x) = [f (x + - f(x)] /b and 

Vmf(x) — \f(x) — / (x — en\a/2\] /b, where b is the dis- 
tance between the nearest neighbors. After some algebra 
calculation, the prefactors A, B, and C are related to the 
usual elastic constant of the crystal /x, A and £ defined 
by the continuum elastic energy 
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and can be expressed for fee as follows: 
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and for bec as 
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Uie l m their components along the link directions. 



For very low temperatures, the atomic positions devi- 
ate very little from those of an ideal crystal. It is there- 
fore suggestive to use these small deviations for defining 
the displacement field. But this definition can be con- 
sistent only for a small time span. Due to fluctuations, 
thermals as well as quantum, the atoms are capable of 
exchanging positions with their neighbors and migrate, 
after a sufficiently long time, through the entire crystal. 
This process of self-diffusion makes it impossible to spec- 
ify the displacement field uniquely. Thus, as a matter of 
principle, the displacement field is undetermined up to an 
arbitrary lattice vector, it is impossible to say whether 
an atom is displaced by un\ (x) or by 

where Nm.(x) is the jumping field with integer value. 
Because of this, we should introduce an extra sum over 
an integer- valued field ng, m ) (x) with a gauge fixing term 
<fr[n(2 m )] in the corresponding partition function, thus the 
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partition function of the defect melting model reads 
dui (x) 
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The partition function is invariant under the following 
defect gauge transformations 

U(i) -> um(x) + bNm(x), n {im ) -> n (im) + WmN (m )(x). 

This partition function can be expressed in terms of the 
stress field <7jj by 
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e \l) e \m) ai i ^ 0r ^ CC lattice, and (3 — 
(a 3 /x)/(nfcsT(27r) 2 ). In the calculation process from (JJJ 
to (jHJ, the corresponding expressions of A, B, C and 
are used. In spite of the tedious intermediate calculation, 
the prefactors in the partition functions JSJ) are the same 
for all cubic lattices. 

For the-lowest order approximation, we need only to 
investigate the system in low-temperature and high- 
temperature limits. It has been shown in the textbook 
(see Fig. 12.1 on p. 1084) that the intersection of 
the free energies of the two limiting curves yields good 
estimates for the melting temperatures in simple-cubic 
lattices. 

For T — > 0, the defect configuration is completely 
frozen out, and the partition function has a classical 
limit. From the lattice energy we obtain the clas- 
sical partition function of cubic lattice 
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where the dimensionless parameter t is defined by the 
trace log 
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and Km and Km are lattice momenta, for fee 

V2 
ai 

and for bec 
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Thus we find the free energy density of the fee and bec 
lattices in the low-temperature limit 
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In the opposite limit of high temperature, the defects 
are prolific and the sum over nn m ) can be approximated 
by integral enforcing (Ty = 0. Then the partition function 
© yields the free energy density 
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From the intersection of this with the low-temperature 
approximation (|13f) . we obtain the melting relation for 
the fee and bec lattices 
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The formula is structurally very similar to the Linde- 
mann rule L = const., but predicts, in addition, the size 
of L to be averagely of the order of 126. Moreover, from 
our expression, one can recognize clearly that the melting 
temperature will be zero when the anisotropic parameter 
£ equals to zero, this is in contrast to recent phenomeno- 
logical theories of melting by other authors [j^. The 
novelty of this formula with respect to a similar formula 
in Ref. 0] [see Eq. (12.13) on p. 1079] lies in a determi- 
nation of trace log term I for bec and fee lattices. The 
resulting melting temperatures from Eq. Ijl5|l are shown 
in Table 1, where they are compared with experimental 
numbers and with the results found from Lindemann's 
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rule where the Lindemann number is initially undeter- 
mined and must be extracted from the average over all 
materials. Our lowest-order theory has an average preci- 
sion of about 12%, which is better by a factor 2 than the 
22% of the numbers from derived Lindcmann's rule. 



Element 


M 


A 




^m,thcor . 


^m,cxp. 


k(%) 


S(%) 


Ag (fee) 


51.10 


97.30 


0.335 


1246.5 


1234.0 


1 


24 


Au (fee) 


45.40 


169.70 


0.351 


1201.4 


1336.2 


-10 


6 


Ba (bec) 


9.50 


8.00 


0.242 


875.4 


998.0 


-12 


22 


Ca (fee) 


16.30 


18.20 


0.294 


889.0 


1112.0 


-20 


-14 


Co (fee) 


128.00 


160.00 


0.320 


1892.7 


1765.0 


7 


-2 


Cu (fee) 


81.80 


124.90 


0.314 


1299.8 


1356.0 


-4 


16 


Li (bee) 


10.80 


12.50 


0.110 


450.2 


454.0 


-1 


-14 


Nb (bee) 


28.70 


134.00 


1.951 


2778.2 


2741.0 


1 


-39 


Ni (fee) 


124.70 


147.30 


0.398 


2034.9 


1726.0 


17 


16 


Pb (fee) 


14.40 


39.20 


0.222 


507.4 


600.6 


-15 


26 


Pd (fee) 


71.20 


176.10 


0.407 


1723.0 


1825.0 


-6 


7 


Pt (fee) 


76.50 


250.70 


0.627 


2527.3 


2042.0 


23 


29 


Sr (fee) 


9.90 


10.30 


0.253 


868.0 


1045.0 


-17 


11 


Ta (bee) 


81.80 


157.40 


0.627 


3879.0 


3271.0 


18 


-32 


Th (fee) 


47.80 


48.90 


0.278 


1882.8 


2024.0 


-7 


43 


Tl (bee) 


11.00 


34.00 


0.309 


609.9 


576.0 


5 


8 


V (bee) 


46.00 


119.40 


1.228 


2444.3 


2178.0 


12 


-15 


W (bee) 


163.10 


204.90 


1.005 


3497.4 


3653.0 


-4 


3 



Table 1. Theoretical results of melting temperatures 
(units in K) derived from the elastic constants of cu- 
bic crystals (units in GPa, gotten from 31]) compared 
with the experimental data (units in K, gotten from 
H3)- The second-last column shows the relative error 
K = (T^thcor/TVcxp - 1) • 100 of our theory which 
is by a factor 2 lower than the relative error 5 = 



m,Lind/ 



/7m iexp — 1) • 100 found from Lindcmann's rule. 



The Lindemann melting temperatures T mj Lind are from 
Ref. El. 



There is no problem, in principle, to calculate system- 
atic improvements to formula (|15|l by including defect 
excitations into the low-temperature approximation l|13(l 
and stress corrections into the high-temperature approxi- 
mation (|14f> . This was done in Ref. Q for the unphysical 
simple cubic lattices, and will be done for bec and fee 
lattices in future work. 
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